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Abstract 

We have studied the multipole polarizability of a graded spherical particle in 
a nonuniform electric field, in which the conductivity can vary radially inside 
the particle. The main objective of this work is to access the effects of mul- 
tipole interactions at small interparticle separations, which can be important 
in non-dilute suspensions of functionally graded materials. The nonuniform 
electric field arises either from that applied on the particle or from the local 
field of all other particles. We developed a differential effective multipole mo- 
ment approximation (DEMMA) to compute the multipole moment of a graded 
spherical particle in a nonuniform external field. Moreover, we compare the 
DEMMA results with the exact results of the power-law graded profile and 
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the agreement is excellent. The extension to anisotropic DEMMA will be 
studied in an Appendix. 

PACS Number(s): 77.22.-d, 77.84.Lf, 42.79.Ry, 41.20.Cv 
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I. INTRODUCTION 



In functionally graded materials (FGM), the materials or microstructure properties may 
vary spatially to meet the specific needs in various engineering applications [1] . The material 
or microstructure gradients in FGM lead to properties quite distinct from those of the 
homogeneous materials and conventional composite materials [1,2]. A unique advantage of 
using FGM is that one can tailor the materials properties via the design of the gradients. 
Over the past few years, there have been overwhelming attempts, both theoretical and 
experimental, of studying the responses of FGM to mechanical, thermal, and electric loads 
and for different microstructure in various systems [1-8]. Such FGM may exist in Nature 
or they can be made artificially. For instance, graded morphogen profiles can exist in a 
cell layer [4]. Graded structure may be produced by using various approaches, such as a 
three-dimensional X-ray microscopy technique [7] , deformation under large sliding loads [8] , 
and adsorbate-substrate atomic exchange during growth [6]. Moreover, it has been reported 
recently that the control of a compatibility factor can facilitate the engineering of FGM [5] . 

There have been various theoretical attempts to treat the composite materials of homo- 
geneous inclusions [9] as well as multi-shell inclusions [10-13]. These established theories 
for homogeneous inclusions, however, cannot be applied to graded inclusions directly. To 
this end, we have recently developed a first-principles approach for calculating the effective 
response of dilute composites of graded cylindrical inclusions [14] as well as graded spherical 
particles [15]. More recently, the spectral representation of graded composites has also been 
established [16]. These theories essentially assumed that the particles are sufficiently far 
apart so that it is possible to neglect contributions from higher-order multipoles and these 
approaches are only valid for dilute composites. As the particles become closer, the local 
field is extremely inhomogeneous on the surface of the particles. To this end, the electro- 
magnetic response of a system of homogeneous (non-graded) spherical particles has been 
studied extensively in the literature [17,18]. In this work, we aim at studying the multipole 
polarizability of a graded spherical particle in a nonuniform field which occurs naturally 
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for a system of interacting particles. The electrostatic boundary-value problem of a graded 
spherical particle will be solved to obtain exact analytic results for the power-law graded 
profile. For arbitrary graded profiles, we will develop a differential effective multipole mo- 
ment approximation (DEMMA) to compute the multipole moment of a graded spherical 
particle to capture the multipole response of the particle in a nonuniform field. 

The plan of the paper is as follows. In section II, we will solve the boundary-value 
problem of a graded spherical particle in a nonuniform electric field. In this way, the exact 
analytic expression for the multipole polarizability is obtained for the case of a power-law 
conductivity profile. In section III, we derive the DEMMA for the multipole polarizability 
of a graded spherical particle of an arbitrary graded profile. In section IV, we compare 
the DEMMA results with the exact results of first-principles approach. Discussion of our 
approaches will be given in section V. The generalization to anisotropic DEMMA will be 
discussed. 

II. FIRST-PRINCIPLES APPROACH 

We consider a graded spherical inclusion of radius a subjected to a nonuniform electric 
field of a point charge placed on the z-axis. For electrical conductivity, the constitutive 
relations read J = Ui{r)E and J = a m E respectively in the graded spherical inclusion and 
the host medium, where <7j(r) is the conductivity profile of the graded spherical inclusion 
and a m is the conductivity of the host medium. Moreover, 

V-J = 0, Vx£ = 0. 

To this end, E can be written as the gradient of a scalar potential $, E — — V$, leading to 
a partial differential equation: 

V • [cr(r)V$] = 0. (1) 

where <r(r) is the dimensionless conductivity profile, while <r(r) = o"j(r)/<7 m in the inclusion, 
and cr(r) = 1 in the host medium. Without loss of generality, we may also let a = 1. 
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In spherical coordinates, the electric potential $ satisfies 

We place a point charge of magnitude g on the z-axis at a distance i? from the center of 
the graded sphere. Thus $ is independent of the azimuthal angle (p. If we write <3> = f(r)Q(9) 
to achieve separation of variables, we obtain two distinct ordinary differential equations. For 
the radial function f(r), 

|:(^)|)-/(/ + iMr)/ = o, (3) 

where / is an integer, while 0(0) satisfies the Legendre equation [9]. The potential can be 
obtained by solving Eq.(3). Exact analytic results can be obtained for a power-law profile 
[14,15], linear profile [14], and exponential profile [19-21]. 

Let us consider the case in which the conductivity profile of the particle has a power-law 
dependence on the radius, a(r) = cr k , with k > where < r < I. Then the radial equation 
becomes 

+ k + 2dl _ /(/ + !)/ = Q 

dr 2 r dr r 2 

As Eq.(4) is a homogeneous equation, it admits a power-law solution [15], 

f(r) = r s . (5) 

Substituting it into Eq.(4), we obtain the equation s 2 + s(k + 1) — 1(1 + 1) =0 and the 
solution is 



k 1 
S =2 



(6) 



-(fc + 1) ± V(/c + l) 2 + 4Z(Z + 1) 
The potentials in the inclusion and the host medium are given, respectively, by 

oo 

<5> l (r,9) = Y,Ar s + (l) Pi(cos9), 

1=0 

n °° 

$ m (r, 0) = - + E Ar- (m) PKcos 9). (7) 



1 1=0 



The positive root of Eq.(6) has been chosen for nonsingular solution inside the inclusion. 
Meanwhile, the potential functions should satisfy the boundary conditions, as follow: 

$i(r,6) \ r =l = $m(r,e) \ r =l , 

°(r)^—U =— ^U- (8) 

The potential due to a point charge located at a distance R from the center of the particle 
can be rewritten in a multipole expansion in the vicinity of the particle surface [9] : 

f = ?£^(cos0). (9) 
' q 1=0 JX 

From Eq.(8), we obtain a set of simultaneous linear equations 

q rj 4 _ ( q l + 1 fl \ 



1 R l + l J ' 1 ^(1)4(0 Vi? m i 
Solving the above equations, we obtain the coefficients: 

q 21 + 1 



A = 



R l + l l(F t + 1) + 1 



q 1{F\ - 1) 

~~ ~lV+^l(Fi + 1) + 1 ' l ~ h (10) 



where 



<j m I 

The coefficients B\ are just proportional to the multipole response of the graded spherical 
particle to the applied multipole field. We thus identify the multipole factor 

H l = 1 KFl '^ ) , 1>1. (12) 
Z(F, + 1) + 1' " V ; 

From Eq.(ll) and Eq.(12), the quantity Fi can be identified as the /-dependent equivalent 

conductivity of the graded particle [see Eq.(16) in section III below]. In the case of a 

homogeneous sphere, k — 0, s+(l) — I, F\ — <7i/<J m , one recovers the well known result [17]: 

l(ai - a m ) 
m = 77 — ■ — , / > l. 

For a uniform field, however, the dipole factor of Dong et al [15] recovers. 

6 



III. DIFFERENTIAL EFFECTIVE MULTIPOLE MOMENT APPROXIMATION 



In this section, we develop the differential effective multipole moment approximation 
(DEMMA) for a graded spherical particle and hence compute the multipole factor. To 
establish the differential effective multipole moment theory, we mimic the graded profile by 
a multi-shell construction [22], i.e., we build up the conductivity profile gradually by adding 
shells. We start with an infinitesimal spherical core of conductivity o"j(0 + ) and keep on 
adding spherical shells of conductivity given by <7j(r) at radius r, until r = a is reached. 
At radius r, we have an inhomogeneous sphere with certain multipole moment. We further 
replace the inhomogeneous sphere by a homogeneous sphere of the same multipole moment 
and the graded profile is replaced by an effective conductivity (Xj(r). Thus, 

#«(r) = „?; (r) ~r m) • (13) 

Next, we add to the sphere a spherical shell of infinitesimal thickness dr, of conductivity 
<7j(r). The resulting multipole factor H[ will change according to [18] 

1 - (Tm/aiirfia^rWaiir) +la t (r)) + pai(r)[l' (a^r) - a^r)) + la m (&i(r) / a^r) - 1)] 



[1 + / / o- m /(/o- i (r))](Ti(r)(/ / (Ti(r) + Icr^r)) + l'p<Ti(r)(cr m - cr^r) + a^r) - (J m <Ji{r) / <Ji{r)) 

(14) 



with I' = 1 + 1 and p = [r/(r + dr)] 2l+1 . 

Of course, the equivalent conductivity &i(r), being related to Hi(r), should also change 
by the same token. Let us write H[ = H[ + dH[, and take the limit dr — > 0, we obtain a 
differential equation: 

^ - - {2l+1 y ffma{r) W{r) + / + H l{ r)l)a m + (H l{ r) - l M r)] 

x [(Htir) + I + mir)^ - (H^r) - l)(l + l)a(r)}. (15) 

Thus the multipole factor of a graded spherical particle can be calculated by solving the 
above differential equation with a given graded profile <7j(r). Solving Hi(r) from Eq.(15) gives 
numerical results for the multipole factor. The nonlinear first-order differential equation can 



be integrated if we are given the graded profile <7j(r) and the initial condition Hi(r = 0). In 
general, it works for arbitrary graded profiles <7j(r). 

The substitution of the relation Eq.(13) into Eq. (15) yields the differential equation for 
the equivalent conductivity <7j(r), 

<&i(r) = [^(r) - a t (r)][(l + l)a t (r) + la^r)] 
dr rai(r) 

The / = 1 limit of Eq. (16) is a special case of the Tartar formula, derived for assemblages 
of spheres with varying radial and tangential conductivities [23]. 



IV. NUMERICAL RESULTS 

Figure 1 displays the comparison between the DEMMA [Eq. (15)] and the first-principles 
approach [Eq. (12)], for the model power-law graded profiles <7j(r) = cr k . Here Hi is plotted 
as a function of c, for different k. As c increases, the multipole factor Hi is caused to increase 
monotonically, for different multipole order I. The effect of I can also increase the Hi slightly. 
Nevertheless, increasing k causes the H t to decrease. Interestingly, the excellent agreement 
between the two methods have been shown in this figure. In fact, the DEMMA is valid for 
arbitrary gradation profiles, as mentioned above. However, exact analytic results obtained 
from the first-principles approach often lack except for a few specific graded profiles like 
the power-law profiles. In view of the good comparison shown in this figure, one can safely 
use the DEMMA to treat other graded profiles which might be not possible, or difficult, to 
be solved by using the first-principles approach. To this end, the good agreement between 
DEMMA results and the exact results can be understood by the fact that Eq.(ll) indeed 
solves Eq.(16) for the power-law profile. Moreover, the approximate DEMMA approach has 
been shown to be exact for spherical particles [24]. 

In Fig. 2, we show the equivalent conductivity <7j(r = a) of a graded spherical particle, 
according to the DEMMA [Eq. (16)], for the same model graded profile <Ji(r) = cr k . Eq. (16) 
shows that <7j(r = a) depends only on the multipole order / and the conductivity gradation 
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profile <7j(r). Here a denotes the radius of the particle, which has been set to be unity 
throughout the paper. Similarly, increasing c or I causes the equivalent conductivity to 
increase. However, opposite behavior is obtained for increasing k. 

In fact, the parameter k in the graded profile <7i(r) = cr k measures the degree of inho- 
mogeneity in the graded particle. From Figs. 1 and 2, it is concluded that the presence of 
inhomogeneity in the particle can affect the electrical properties of the particles significantly. 
In other words, once an inhomogeneous particle is used in reality, its effect of inhomogeneity 
should be taken into account. 

V. DISCUSSION AND CONCLUSION 

Here a few comments are in order. In this work, we have developed a differential effective 
multipole moment approximation (DEMMA) to compute the multipole moment of a graded 
spherical particle. We compared the DEMMA results with the exact results of the power-law 
profile and the agreement is excellent. Note that an exact solution is very few in composite 
research and to have one yields much insight. Such solutions should be useful as benchmarks. 

We are now in a position to propose some applications of the present theory. As the 
multipole response is sensitive to the graded profile of the particles as well as to the structure 
of the nonuniform field source, there is a potential application to ac electrokinetics of graded 
colloidal particles [25]. The similar approach can be applied to electrorheological (ER) fluids 
[26,27] because the particles in ER fluids can have very dense structures locally. The local 
electric fields are extremely inhomogeneous near the particles so that multipole effects can 
play an important role. In this regard, we can also study the interparticle force between 
graded particles [28] in ER fluids. 

In the other topics, we may attempt the similar calculation of the multipole response 
of a graded metallic sphere in the nonuniform field of an oscillating point dipole at optical 
frequency. The graded Drude dielectric function can be adopted [29]. When the oscillating 
source is placed close enough to the graded metallic sphere, higher-order multipole response 
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can be excited. The approach may also be applied to the electroencephalogram of the human 
brain by regarding the brain as a graded anisotropic conducting sphere [30]. In this way, 
the forward problem can be analyzed. The derivation of the anisotropic DEMMA will be 
given in Appendix A. 

In summary, we have studied the multipole response of a graded spherical particle in a 
nonuniform field. We developed a DEMMA to compute the multipole moment of a graded 
spherical particle. Moreover, we compare the DEMMA results with the exact results of the 
power-law profile and the agreement is excellent. 
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APPENDIX A: ANISOTROPIC DIFFERENTIAL EFFECTIVE MULTIPOLE 

MOMENT APPROXIMATION 



Let us consider a graded spherical particle with radius a. We adopt the spherical coor- 
dinates for convenience. The graded spherical particle has a tangential conductivity in the 
plane orthogonal to the radial vector of the sphere <J^{r) and a radial conductivity cr^(r). 
Both 0^-(r) and crj'(r) will be prescribed by radial functions, < r < a. In view of the 
symmetry, the anisotropic conductivity of the graded sphere can be expressed as a tensor 
form <jj (r), namely, 

°i (r) = a} (r)ff + a±(r)66 + ^(r)^. (Al) 

Next, we present an anisotropic differential effective multipole moment approximation 
(ADEMMA), which, similar to DEMMA, is a numerical method for the analysis of the 
electric property of anisotropic graded particles with arbitrary gradation profiles. Similarly, 
we may regard the gradation profile as a multishell construction. In detail, we establish 
the electric profile gradually by adding shells. Let us start with an infinitesimal isotropic 
spherical core with conductivity <Xj(0 + ), and keep on adding shells with both tangential and 
normal conductivity profiles &^(r) and <7j"(r) at radius r, until r = a is reached. At radius 
r, we have an inhomogeneous particle, and further regard such an inhomogeneous particle as 
an effective homogeneous one with an effective conductivity <7j(r), which has the multipole 
factor 

#«(r) = 'f^^ • (A2) 

Then, we add to the particle a shell with infinitesimal thickness dr, with conductivities 
o" i - L (r) and <7j"(r). Its multipole factor H[ should change according to the multipole factor 
of a single-coated particle [31] 

ffl = (o r J l ( r ) M - ~ °i{r)l){(j j \\r)u + - (T rn l) - pi{a}{r)u+ - ai(r)l)(a^(r)u- - a m l) 
(ffi(r)Z - a i W{r)u + ){a i \\{r)u- + a m V)pi - (^i(r)Z - a i W{r)u-)(a i W{r)u + + a m V) ' 
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with u± = [-1 ± yfl + 4/(1 + l)a^(r)/a, i \\(r)]/2, Pl = [r/(r + dr)] u +~ u - , and I' = 1 + 1. Let 
us write further H[ = Hi + dH[, and take the limit dr — > 0, we obtain a differential equation 

^ = - (l + 20rt mff< ll(r) [fffm " + H ^ l ' a - + u -^"( r ))] 

[Z<r m - w+a^r) + H l {r-){l'a m + u+a^r))]. (A4) 

Thus the multipole factor of an anisotropic graded spherical particle ifj(r = a) can be 
calculated by solving the nonlinear first-order differential equation [Eq. (A4)] which can be 
integrated, at least numerically if we are given the graded profiles cr^ir) and a,^(r) and the 
initial condition Hi(r = 0). The substitution of Eq. (A2) into Eq. (A4) yields the differential 
equation for the equivalent conductivity 

da;(r) (1 + /)crJI(r)o- i ± (r) - orjl(r)o-i(r) - /^(r) 2 

~1\ — = • ( A5 > 

dr rai"{r) 

Eqs. (A4) and (A5) can respectively reduce to Eqs. (19) and (20), as long as there is o^(r) = 
<7;"(r) = <7j(r). The substitution of I — 1 into Eq. (A5) reduces to the Tartar formula, derived 
for assemblages of spheres with varying radial and tangential conductivities [23]. Last but 
not least, it is also instructive to extend the first-principles approach of Ref. [20] to multipole 
response and compare with the ADEMMA results. 
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FIGURES 

FIG. 1. Comparison of multiple factor Hi, between the DEMMA [Eq. (15)] and the 
first-principles approach [Eq. (12)]. The graded profile is <7j(r) = cr k and a m = 1. Excellent 
agreement between the two methods have been shown. 

FIG. 2. The equivalent conductivity <7j(r = a) of a graded spherical particle, according to the 
DEMMA [Eq. (16)]. Here a denotes the radius of the particle, which is set to unity throughout 
the paper. The graded profile is <7j(r) = cr k . 
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